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Abstract 


Controlled variation of the electronic properties of 2D materials by applying strain has 
emerged as a promising way to design materials for customized applications. Using first prin¬ 
ciples density functional theory calculations, we show that while the electronic structure and 
indirect band gap of SnS 2 do not change significantly with the number of layers, they can be 
reversibly tuned by applying biaxial tensile (BT), biaxial compressive (BC), and normal com¬ 
pressive (NC) strains. Mono to multilayered SnS 2 exhibit a reversible semiconductor to metal 
transition (S-M) at strain values of 0.17, —0.26, and —0.24 under BT, BC, and NC strains, 
respectively. Due to weaker interlayer coupling, the critical strain value required to achieve 
S-M transition in SnS 2 under NC strain is much higher than for M 0 S 2 . The S-M transition for 
BT, BC, and NC strains is caused by the interaction between the S-p^ and Sn-5, S-px/py and 
Sn-5, and S-p^ and Sn-5 orbitals, respectively. 


Introduction 

The advent of semiconducting two-dimensional ( 2 D) materials such as, MoS2,-”^ phosphorene,- 
and SnS2-'"- has attracted immense interest in the scientific community because of their poten¬ 
tial applications in electronic, magnetic, optoelectronic, and sensing devices. Modifying their 
electronic properties in a controlled manner is the key to the success of technologies based on 
these materials. Some popular techniques used in modifying their properties are electric field, 
heterostructuring, functionalization, and application of strain. Due to experimental feasibility of 
strain, it has been a very effective way of transforming and tuning electronic mechanical— 
as well as magnetic— properties of these materials. Recent studies have shown its applicabil¬ 
ity in modifying some of the fundamental properties of bilayer phosphorene and transition metal 
dichalcogenides (TMDs) where a complete semiconductor to metal (S-M) transition was observed 
under the applied strain.-"— The effect of strain has been studied in not only two-dimensional ma¬ 
terials, but also, in materials of other dimensionalities. Some examples include direct to indirect 
band gap transition in multilayered WSe2— and drastic changes in the electronic properties of few 
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layers to bulk M0S2.— Its application is not only limited to electronic properties but it can also 
engineer the vibrational properties-^*^!^ of these materials. 

Among all the layered 2 D materials, semiconducting SnS2 is more earth-abundant and environment- 
friendly and demonstrates a wide spectrum of applications. For example, lithium-^-— and sodium- 
ion batteries— with anodes constructed from SnS2 show high capacity, enhanced cyclability, and 
excellent reversibility. In addition, SnS2 shows good photocatalytic activityNanosheets made 
from SnS2 are considered to be very good for hydrogen generation by photocatalytic water split- 
ting.-i 22 The visible-light photocatalytic activity of reduced graphene oxide has been shown to be 
enhanced by doping Cu in SnS2 sheets.— Further usage of SnS2 in photonics has been demon¬ 
strated through fabrication of fast-response SnS2 photodetectors having an order of magnitude 
faster photocurrent response times— than those reported for other layered materials. This makes 
SnS2 a suitable candidate for sustainable and “green" optoelectronics applications. However, most 
of these applications rely on the inherent band gap of SnS2 and therefore, its functionality could 
be further enhanced by tuning its band gap. Here, we explore this possibility for layered SnS2 as 
a function of (i) number of layers and (ii) application of different strains such as biaxial tensile 
(BT), biaxial compressive (BC), and normal compressive (NC) strains. Band gap tuning in lay¬ 
ered materials such as 14082—*^ and phosphorene^*^ has been achieved using these strategies. 
Interestingly, at zero strain, on increasing the layer number in SnS2, the band gap was found to be 
indirect and insensitive to the number of the layers due to the weaker interlayer coupling in SnS2 
compared to M0S2. In addition, applying strain does not change the nature of the band gap. On 
the contrary, irrespective of the type of applied strain, a reversible semiconductor to metal (S-M) 
transition was observed. The S-M transition in bilayer ( 2 L) SnS2 was achieved at strain values of 
0 . 17 , — 0 . 26 , and — 0 . 24 , for BT, BC, and NC strains, respectively. The strain values required to 
achieve this S-M transition in SnS2 are higher in magnitude compared to M0S2, which is attributed 
to the weak interlayer coupling in SnS2. 
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Methodology 


The calculations were performed using first-principles density functional theory (DFT)^ as im¬ 
plemented in the Vienna ab initio simulation package (VASP).— Projector augmented wave 
(PAW )^>22 pseudopotentials were used to represent the electron-ion interactions. In order to 
obtain an accurate band gap, the Heyd-Scuseria-Ernzerhof (HSE 06 ) hybrid functional—”^ was 
used, which models the short-range exchange energy of the electrons by fractions of Eock and 
Perdew-Burke-Ernzerhof (PBE) exchange.— The addition of the Eock exchange partially removes 
the PBE self-interaction which resolves the band gap overestimation problem. Sufficient vac¬ 
uum was added along the z-axis to avoid spurious interactions between the sheet and its periodic 
images. To model the van der Waals interactions between the SnS2 layers, we incorporated a 
semi-empirical dispersion potential (D) to the conventional Kohn-Sham DET energy, through a 
pair-wise force field through Grimme’s DET-D 3 method.— The Brillouin zone was sampled by a 
well-converged 9 x 9 x 1 Monkhorst-Pack— k-grid. Structural relaxation was performed using the 
conjugate-gradient method until the absolute values of the components of the Hellman-Eeynman 
forces were converged to within 0.005 eV/A. 

Results and discussion 

SnS2 belongs to the group-IV dichalcogenide MX2 family,— where, M and X are group IV and 
chalcogen elements, respectively. The hexagonal SnS2 monolayer has a Cdl2 layered crystal struc¬ 
ture with a space group symmetry of P 3 ml— unlike M0S2, which is a transition metal dichalco¬ 
genide with a space group symmetry of P 6 T,/mmc. It consists of three atomic sublayers with 
Sn atoms in the middle layer covalently bonded to S atoms located at the top and bottom layers 
(Pig. [Ha)). In multilayered SnS2, the interlayer interactions are van der Waals type.— “— The op¬ 
timized lattice parameters of bulk and monolayer SnS2 are abuik = 3.689 A, Cbuik = 5.98 A, and 
«mono = 3.69 A, respectively, which are in good agreement with the experimental values.-^^^ SnS2 
prefers AA stacking in the multilayered structure (Pig. Hb)) and the in-plane lattice parameter a 
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and interlayer distanees remain the same as in bulk. 

Bulk SnS 2 is an indirect band gap semiconductor with an HSE06 band gap of 2.27 eV. The 
valence band maxima (VBM) is in between the F and M points, while the conduction band minima 
(CBM) is at the L point. The band structures of 1, 2, 3, 4, and 5 layered SnS 2 are calculated and 
are shown for 1-3 layers in Fig. [IJc). The band structures of all these multilayers show that the 
number of bands forming VBM and CBM increase with layer number and are equal to the number 
of layers, as shown in Fig. [He). This gives an indication of weak interlayer interaction. Fike 
bulk, all these multilayered SnS 2 structures are indirect band gap semiconductors. The HSE06 
band gap values of 1, 2, 3, 4, and 5 layers are 2.47 eV, 2.46 eV, 2.40 eV, 2.35 eV, and 2.34 eV, 
respectively. The overall variation in the band gap from monolayer to bulk is very small (0.20 eV) 
demonstrating insensitivity of the electronic structure to the number of layers. It is well-known that 
experimental growth of these layered materials with a precise control over (a) the number of layers 
and (b) the quantity of material of a pre-decided thickness, is a difficult task. Therefore, although 
this insensitivity will limit its application in optical and sensing devices, it can still be beneficial 
from the nanoelectronics perspective. The devices fabricated from SnS 2 may have minimum noise 
or error in terms of the change in (a) the contact resistance and (b) carrier mobilities, both arising 
from the difference in the number of layers, unlike for M 0 S 2 field effect transistors (FETs).— 
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Figure 1: Structures of unstrained (a) monolayer SnS 2 and (b) multilayered SnS 2 of equilibrium 
layer thickness Jq- (c) Band structure of unstrained one- (IF), two- (2F) and three-layer (3F) SnS 2 . 
Fermi level is set to 0. 

Having shown the insensitivity of the band gap towards the number of layers, we now explore 
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the possibility of electronic structure tuning by applying strain. For any strain value, the nature 
of the band gap (i.e. indirect band gap) was found to be independent of the number of layers. 
Therefore, here, we will discuss only the results of bilayer (2L) SnS 2 in detail. First, we study 
the effect of in-plane uniform BC and BT strains on the electronic structure of 2L-SnS2. The 
in-plane strain (e) is given by e = where, a and ao are lattice parameters with and without 
strain, respectively. For the unstrained 2L-SnS2, the VBM is in between the F and M points, 
while the CBM is at the M point. The conduction bands CBM and CBM—1 are dispersive at 
the M point, whereas, the valence bands VBM and VBM—1, are not very dispersive. Hence, the 
overall curvature as well as the mobilities of these valence bands are expected to be quite low. 
With increasing BT strain, the VBM and VBM—1 remain flat and are still located in the F-M 
region as shown in Fig. [2ta). At a strain value of 0.06, the CBM shifts to the F point while still 
preserving the indirect nature of the band gap. In addition, the curvature of the conduction bands at 
the M (F) point decreases (increases), implying that, locally, the mobility of n-type carriers at this 
point may decrease (increase). However, the overall mobility of the n-type carriers may not change 
significantly, owing to the compensation of these local changes in the mobilities. Also, irrespective 
of the strain, the shape of the conduction as well as the valence bands do not change. At a critical 
strain value of 0.17, the CBM at F crosses the Fermi level (Fig. [If a)), rendering 2L-SnS2 metallic 
(Fig.iab)). 

Under BC strain, similar changes in the band dispersion of 2L-SnS2 are observed. Surprisingly, 
unlike the BT strain, at a strain value of —0.06, there is a significant shift in the CBM away from 
the Fermi level as compared to the VBM resulting in an increase in the band gap. However, 
irrespective of strain, the CBM remains at the M point. The VBM remains at the F-M location 
upto a strain of —0.06, beyond which, it moves to the F point. In addition, both, the valence as 
well as conduction bands become more dispersive with strain, indicating that the overall mobilities 
of p- and n-type carriers may increase. Similar to BT strain, the movement of VBM and CBM 
towards the Fermi level leads to the band gap reduction, while still preserving an indirect band 
gap, as shown in Fig.[2tc). At a critical strain of —0.26, the S-M transition occurs (Fig.[2td)), with 
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the VBM crossing the Fermi level (Fig. [2c)). 




Figure 2: Band structures of bilayer SnS 2 as a function of (a) BT, (c) BC, and (e) NC strains. Fermi 
level is set to 0. Variation of band gap from IL-to 5L-SnS2 as a function of (b) BT, (d) BC, and (f) 
NC strains. 

We next study the changes in band structure of SnS 2 under normal compressive (NC) strain 
(Cnc) defined as £nc = where d and do (Fig. Etb)) are the layer thickness at applied and 
zero strains. The minimum required number of layers for NC strain is two. The NC strain was 
applied perpendicular to the plane of the multi-layers of SnS 2 . Constrained relaxation mechanism 
was incorporated wherein the atoms in upper and lower layers were restricted to move along the 
normal direction at each strain value. With increasing strain, the band gap reduces smoothly and 
becomes metallic at —0.24, as shown in Figs.[2e) and (f). In the strain range of —0.08 to —0.17, 
the VBM remains in between the F and M points, whereas there is a drastic shift in the position 
of CBM from M to in between F and M points. Similar to BC strain, the bands become more 
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dispersive with increasing NC strain, indicating that the overall mobilities of p- and n-type carriers 
can increase. At —0.24 strain, a band crossing occurs at the Fermi level rendering it metallic. 

Next, we discuss the mechanism of S-M transition under BT, BC, and NC strains. We cal¬ 
culated the angular momentum projected density of states (LDOS) and band-decomposed charge 
density of 2L-SnS2 as a function of these strains. In unstrained SnS 2 , the valence and conduction 
bands originate predominantly from S-px/py/pz, and the direction-independent Sn-5 orbitals, re¬ 
spectively, as shown in Fig. [2 a). With increasing BT strain up to 0.06, S-pz crosses the S-px and 
S-py orbitals, and becomes the frontier orbital. The band gap begins to reduce due to the move¬ 
ment of these frontier orbitals (Sn-5 and S-pz) towards the Fermi level. For strains beyond 0.06, 
the intralayer interactions between these orbitals begin to increase. This is also confirmed by the 
band-decomposed charge density plots (Fig. [3tb)). At a critical strain of 0.17, there is a strong 
hybridization between these orbitals, as is also observed in the CBM band-decomposed charge 
density (Fig. Oc)), which finally leads to the semiconductor to metal transition. 

In the case of BC strain, the S-pz orbital crosses the tails of the S-px and S-py orbitals and 
moves away from the Fermi level, as shown in Fig. Od). Hence, the S-px/py and Sn-5 orbitals 
contribute predominantly to the valence and the conduction bands, respectively. With a strain upto 
—0.06, the band gap initially increases, as observed from the band structure in Fig. [I^d). This 
increase is due to the slight upward (downward) shift in the Sn-5 {S-pJpy) orbitals along with the 
crossing of the S-pz orbital, causing the conduction (valence) bands to move up (down) as shown 
in Fig. [2d)- A similar band gap increase at an intermediate BC strain is observed in M 0 S 2 .— As 
shown in the VBM band-decomposed charge density in Fig. [2d), the wavefunction spread between 
the S atoms (above and below Sn) in a layer decreases compared to the unstrained case (Fig. [2a)) 
due to the shifting of the valence orbitals px, Py, and Pz of the S atom. This localization causes the 
increase in the band gap. As we go beyond —0.06, the band gap begins to decrease because of the 
movement of the S-pJpy (in-plane orbitals) and Sn-5 orbitals towards the Fermi level. At a critical 
strain of —0.26, this movement results in S-M transition. Here, the in-plane orbitals contribute to 
S-M transition, whereas in BT strain, in addition to the in-plane orbitals, out-of-plane orbitals also 
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Figure 3: (a)-(f) LDOS and corresponding band-decomposed charge densities (VBM and CBM) 
for 2L-SnS2 under BT and BC strains. The isosurface has been set to 0.006 e/A^. 


play a role in triggering S-M transition. 

With increasing NC strain up to —0.17, the Sn-5 and S-pJpylp^ orbitals begin to move towards 
the Fermi level, as shown in Figs. Ufa) and (b). For strains beyond —0.17, these orbitals start over- 
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Figure 4: (a)-(c) LDOS and corresponding band-decomposed charge densities (VBM and CBM) 
for 2L-SnS2 under NC strain. The isosurface has been set to 0.006 dA?. 


lapping because of the reduction in the interlayer distance. At a critical strain of —0.24, interlayer 
interactions become enhanced and as a consequence, results in strong hybridization between S-p^ 
and Sn-5 orbitals (Fig. I5c)). These orbitals cross the Fermi level and lead to the S-M transition. 
From the band-decomposed charge analysis, at this critical strain, the presence of out-of-plane 
lobes (pz) in the VBM and s orbital in CBM are mainly responsible for S-M transition, as shown 
in Fig.llfc). 

Among the applied strains, BT strain results in the fastest S-M transition. Although the critical 
strain required in the case of BT strain is the lowest (|e| = 0.17 compared to |e| = 0.26 and 0.24 in 
the case of BC and NC strains, respectively), the NC strain is more feasible compared to the biaxial 
strains. This is because the intralayer bonding between the atoms is much stronger compared to 
the weaker interlayer van der Waals interactions. 

Having explored band gap tuning using BT, BC, and NC strains, it is interesting to perform a 
comparative study of the band gap variation as a function of (a) interlayer coupling between the 
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Figure 5: (a) Band gap variation as a function of number of layers, (b) interlayer coupling energy 
for the corresponding interlayer distance under NC strain, and (c) band gap as a function of NC 
strain %, for both SnS 2 (red line) and M 0 S 2 (black line). 

layers and (b) applied NC strain, in both 2L-SnS2 and 2 L-M 0 S 2 . Here, we use the generalized 
gradient approximation (GGA) of the PBE form to model the exchange-correlation. Fig. [Sfa) 
shows the band gap as a function of the number of layers for the unstrained structures of M 0 S 2 
and SnS 2 . As we go from IL to 5L, the change in the band gap is small for SnS 2 (~ 0.24 eV) 
compared to that in M 0 S 2 (~ 0.59 eV), indicating that the band gap of SnS 2 is less sensitive to the 
number of layers. The reason for this insensitivity is attributed to the weak interlayer coupling in 
SnS 2 compared to M 0 S 2 . Fig.|3b) shows the interlayer coupling energy as a function of interlayer 
distance for 2L-SnS2 and 2 L-M 0 S 2 . Here, the adjacent layers in SnS 2 are weakly coupled to each 
other, with a coupling energy of 235 meV/unit cell, relative to that of M 0 S 2 (430 meV/unit cell). 
The weaker interlayer coupling in SnS 2 can be explained based on the higher Sn electron affinity 
(Zsn = 107.3 kJ/mol)— compared to Mo electron affinity (;(^mo = 72.3 kJ/mol).— A higher x value 
indicates a lower tendency to interact with the neighboring atoms. Therefore, the polarizability of 
the Sn atom is much lower than in Mo causing a weaker interlayer coupling in SnS 2 . The above 
observation is also confirmed by the charge accumulation and depletion calculations of lL-SnS 2 
and IL-M 0 S 2 as shown in Figs. |6ta)-(d) wherein, the lower accumulation of charge around the 
S atoms and lower depletion of charge around Sn atoms indicates weaker polarizability. As a 
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Figure 6: Charge accumulation (light-green) and depletion (red) of (a)-(b) SnS 2 and (c)-(d) M 0 S 2 . 
The isosurface is set to 0.006 dA?. 

consequence, the out-of-plane linear elastic modulus (C^), i.e. 



( 1 ) 


dz^ 


where, z is the interlayer distance and F’coupiing is the interlayer coupling energy, of SnS 2 (Cz— 11.1 
GPa) is nearly half of that of M 0 S 2 (C^ = 20.6 GPa). As a function of NC strain, in bilayer M 0 S 2 , 
due to stronger interlayer coupling, the strain dependence of band gap is more prominent (0.11 
eV/unit strain) than SnS 2 (0.026 eV/unit strain), as shown in Fig. [S^c). Therefore, the overlap of 
wave functions within the adjacent layers is much stronger in M 0 S 2 than in SnS 2 , and hence a small 
change in the interlayer distance can lead to a band renormalization in M 0 S 2 . For the same reason, 
semiconductor to metal transition is observed at a relatively high strain 23% for 2L-SnS2 compared 
to 2 L-M 0 S 2 (12%) under NC strain (Fig. Sc)). Therefore, the insensitivity of the indirect band 
gap of SnS 2 towards the number of layers and less dependency on strain can be attributed to the 
weak coupling between the layers. 

In conclusion, using density functional theory based calculations, we have shown the tuning 
of the band gap as well as reversibility in S-M transition for multilayers of SnS 2 under different 


12 



strains (BT, BC, and NC). With increasing strain, there is a smooth reduetion in the band gap for 
all the strains, with a small inerease observed at —0.06 BC strain. The eritieal strain at whieh 
S-M transition is aehieved, deereases signifieantly, as we go from mono to five layers for all the 
strains. Moreover, in order to understand the S-M transition under different strains, we analyzed the 
eontribution from different moleeular orbitals by performing LDOS and band-deeomposed eharge 
density ealeulations. For BT, BC, and NC strains, the interaetion between the S-p^ and Sn-5, S- 
Px/Py and Sn-5, and S-p^ and Sn-5 orbitals, respeetively, trigger S-M transition. For eomparison, 
we studied the band gap trends in 2L-SnS2 and 2 L-M 0 S 2 . At zero strain, the ehange in the band 
gap of SnS 2 is ~ 0.24 eV, whieh is more than half of that in M 0 S 2 (~ 0.59 eV), while going 
from monolayer to five layer. Moreover, the band gap in SnS 2 remains indireet irrespeetive of the 
number of layers as well as applied strain. The insensitivity of the band gap in SnS 2 to the number 
of layers is attributed to the weaker interlayer eoupling (235 meV/unit eell) than in M 0 S 2 (430 
meV/unit eell) due to lower polarizability. In terms of band gap dependenee under NC strain, SnS 2 
is less sensitive to strain (0.026 eV per unit strain) than M 0 S 2 (0.11 eV per unit strain). 
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